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Selection of Crystal Chirality: Equilibrium or Nonequilibrium? 
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To study the solution growth of crystals composed of chiral organic molecules, a spin-one 
Ising lattice gas model is proposed. The model turns out to be equivalent to the Blumc- 
Emcry-Grifhths model, which shows an equilibrium chiral symmetry breaking at low tem- 
peratures. The kinetic Monte Carlo simulation of crystal growth, however, demonstrates that 
Ostwald ripening is a very slow process with a characteristic time proportional to the system 
size: The dynamics is nonergodic. It is then argued that by incorporating grinding dynamics, 
homochirality is achieved in a short time, independent of the system size. Grinding limits 
cluster sizes to a certain range independent of system size and at the same time keeps the 
supersaturation so high that population numbers of average-sized clusters grow. If numbers 
of clusters for two types of enantiomers differ by chance, the difference is amplified expo- 
nentially and the system rapidly approaches the homochiral state. Relaxation time to the 
final homochiral state is determined by the average cluster size. We conclude that the system 
should be driven and kept in a nonequilibrium state to achieve homochirality. 

KEYWORDS: homochirality, crystal growth, lattice gas model, grinding, kinetic Monte Carlo 
simulation, driven closed system 



1. Introduction 

Spontaneous symmetry breaking is observed in various fields of physics. In life, the chiral 
symmetry of organic molecules is broken, and the origin of this one-handedness or homochi- 
rality has been a long mystery ever since Pasteur found the molecular chirality. 1 ' 2 The first 
theoretical model for explaining the spontaneous breaking of chiral symmetry in life is pro- 
posed by Frank in terms of a set of rate equations with linear autocatalysis and nonlinear 
mutual inhibition. 3 The actual chemical system that realizes the amplification of enantiomeric 
excess (ee) was found by Soai and coworkers. 4,5 The time development of the yield and ee 
of the products was elucidated by the second-order autocatalytic process, which provides 
nonlinearity. A theoretical model with an additional recycling process was found to achieve 
homochirality. 6 Recently, new chemical systems that show ee amplification have been found. 7 

Chirality also appears in the structure and shape of crystals. In fact, Pasteur found two 
enantiomers of the sodium ammonium salt of tartaric acid, since each enantiomer formed 
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separate crystals with mirror symmetric shapes and different enantiomers were resolved just 
by hand picking. 1 When two enantiomers form separate crystals, the molecule is said to form 
a conglomerate. It is noted, however, that only 10 percent of chiral organic molecules form a 
conglomerate and the rest crystallize into racemic crystals. 

Intrinsic chirality is not a prerequisite for the formation of chiral crystals: Achiral molecules 
can form chiral crystals. The eminent example is quartz, which can be either in the d- or 1- 
form. Sodium chlorate, NaC103, is another example that crystallizes into a chiral form. 8 It 
forms a conglomerate, and d- and 1-crystallites are separately formed. However, by simply 
evaporating a sodium chlorate solution in a growth cell, almost the same amounts of d- and 
1-crystallites are grown, and the original solution ends up into a racemic mixture in the cell. 8 ' 9 
In 2005, Viedma discovered that when the solution is stirred violently and the growing crystal 
is ground to small pieces, all crystals eventually turn into crystals of a single chirality in a 
short time. 10 Later, Noorduin and coworkers 11,12 grew crystals of an organic molecule, the 
imine of 2-methyl-benzaldehyde and phenylglycinamide, under grinding, and they found that 
a molecular chirality converts to a single one. Recently, other organic substances as amino 
acid are found to show similar deracemization phenomena. 13 ' 14 

There have been many theoretical studies on the selection of crystal chirality mainly based 
on rate equations. 15-17 In our recent study, 18 we proposed a simple lattice gas model and 
performed Monte Carlo simulation to compare results with those of Noorduin's experiment. 
Our simulation study revealed that the spontaneous chiral symmetry breaking occurs in a 
short time when one includes grinding and the nonlinear amplification of chirality conversion 
on the surface of enantiomeric crystals. 18 In the present study, we analyze the equilibrium 
and nonequilibrium aspects of the model in detail, and discuss what are the essential features 
of the chiral symmetry breaking in crystal growth. 

In the next section, a lattice model for chiral crystal growth, which has three possible 
states on a lattice site, is presented. The model turns out to be equivalent to the Blume- 
Emery-Griffiths model 19 with chiral symmetry breaking in equilibrium. In §3, we describe 
kinetic processes consisting of crystal growth and dissolution, chirality conversion, and grind- 
ing. Simulation results are described in §4. Without grinding, chirality is selected via the 
Ostwald ripening 20 of the two (or a few) largest clusters of different enantiomers, and it takes 
an enormously long time to establish homochirality with a relaxation time of the order of the 
system size: The system is nonergodic. By including the grinding procedure, the relaxation 
time is shown to be independent of the system size, and the chirality is selected rapidly even 
in a large system. Discussions are given in §5. 

2. Equilibrium Selection of Chirality 

We consider the crystal growth of organic molecules from a solution. The organic molecules 
are chiral and can be in one of the two possible enantiomeric forms, classified as either R or S. 
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It is experimentally known that the grown crystals can show chiral symmetry breaking if they 
form a conglomerate, namely, the R crystal grows separately from the S crystal. Otherwise, 
if chiral molecules form a racemic compound, both R and S enantiomers are contained in 
a single crystallite and achieving homochirality is impossible. This fact indicates that for 
the formation of a homochiral crystal, bonding energies between a pair of the same kind of 
enantiomers, Jrr and J55, should be larger than that between a heteropair, Jrs- 

Jrr = Jss > Jrs- (2-1) 

Due to the microscopic chiral symmetry, bonding energies between the same kind of enan- 
tiomers are assumed to be independent of the enantiomeric type. 

We construct a lattice gas model such that on each lattice site there is either an R or 
S chiral enantiomer, or none. The double occupancy of a lattice site is forbidden. Because 
of simplicity, we consider a two-dimensional system where the space is divided into a square 
lattice. The state of the system changes as molecules change their sites or chirality. When a 
pair of nearest-neighboring sites are occupied by molecules, they form a bond with an energy 
gain. In this section, we consider the equilibrium state of the system at a given temperature 
T. It is expected that at sufficiently low temperatures, enantiomers of the same type will form 
clusters to gain energy, and the chiral symmetry of the system may break spontaneously. 

The model, in fact, reduces to a spin-one Ising model studied by Blume, Emery and 
Griffiths, 19 which is called the BEG model. The state on the site i is described by the spin 
variable Si which, in our model, is assumed to take the value +1 if the site is occupied by an 
R enantiomer, —1 if the site is occupied by an S enantiomer, and if the site is vacant; 

+ 1 occupied by R 
Si = \ vacant (2-2) 
— 1 occupied by S. 
The total number of chiral molecules is given by the sum of Sf as 

N = N R + N s = Y,Sl (2-3) 

i 

whereas the population difference of enantiomers is given by the sum of Si as 

M = N R -N s = J2 S i- (2-4) 

i 

Here, Nr and N$ represent the numbers of R and S enantiomers, respectively, and the sum- 
mation runs over a square lattice of L 2 sites. The average value of Sf corresponds to the 
concentration of chiral molecules as 

c=(Sf) = (N)/L 2 , (2.5) 
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and the average value of Si is related to the crystal enantiomeric excess (cee) or the chiral 
order parameter <f> as 

aj> = (Si) = (M)/L 2 . (2.6) 

One can construct projection operators for each state as follows: Pi(R) = Si(Si + l)/2 
is a projection operator for an R enantiomer since it is unity if the site i is occupied by an 
R enantiomer, and zero otherwise. Similarly, Pi(S) = Si(Si — l)/2 and Pj(0) = 1 — Sf are 
projection operators for an S enantiomer and a vacancy, respectively. By the use of these 
projection operators, the interaction energy is written as 

E o = - + J R s(Pi(R)Pj(S) 

(id) 

+ PifflPjiR)) + JssPi(S)Pj(S)] 

= - Y,( KS ? S j + js i s j) ( 2 - 7 ) 

(id) 

with 

K=^(Jrr + Jrs), J =\(Jrr-Jrs)- (2.8) 

Here, the symmetry relation Jrr = J$s is utilized. To discuss phase transitions, we introduce 
conjugate intensive variables, namely, chemical potentials {j,r and /is, as 

E = E Q - urNr - n s N s 

= - ^(KS^I + Js . Sj) -^tf-H^Si (2.9) 

{id) i i 

with \x = (fJ,R + //s)/2 and H = (/ir — /xg)/2. This is just the model discussed by Blume et 
al. 19 The field H breaks the chiral symmetry externally. The spontaneous chiral symmetry 
breaking will take place in the absence of an external field, or H = 0. 

For a macroscopic state defined by the numbers of R and S enantiomers, Nr and Ng, 
respectively, the energy in the mean field approximation is written as 

E = - ^{JrrN 2 r + J ss Nl + 2Jr S NrN s ) 

-fi(N R + N s )-H(N R -Ns) (2.10) 
with the coordination number z = 4 for a square lattice, and the entropy as 

5 = fc B ln -. (2.11) 

N R \N S 1(L 2 - N R - N s )l y ' 

For simplicity, we take the Boltzmann constant to be unity, ks = 1, hereafter. In terms of 
the molecular concentration c = (Nr + Ng)/L 2 and the chiral order parameter <\> = (Nr — 
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By minimizing / in terms of c and (ft & t the fixed intensive parameters ji and H , one obtains 
the self-consistent equations for equilibrium values of c and (ft as 

zJc(ft + H 



tanh ■ 



T 



2 cosh[(z Jc<j> + H)/T] 



" exp[-(zKc + fi)/T\+2coBh[(zJat> + H)/T\' ^ 2 ' 13 ^ 
To study the spontaneous chiral symmetry breaking, the external field H is set to be zero, 
whereas the chemical potential [i is chosen appropriately to fix the concentration c at a desired 
value. 




For H = 0, there are two solutions; one is racemic with <f> = 0, and the other is chiral 
with eft ^ 0. The true equilibrium state is that with the minimum free energy /. For the 
case K = where Jrs = —Jrr < so that R and S are repulsive, the phase diagram is 
shown in Fig. 1(a), which is essentially the same as that depicted by Blume et al. 19 (Note 
that the He 3 concentration x in ref. 19 is related to our molecular concentration c by the 
relation x = 1 — c.) At high temperatures, T > zJ, the system is racemic for all values of 
concentration c. Below the second-order transition temperature zJc, the chiral symmetry is 
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broken with a finite value of ^ at a given c. Below the tricritical temperature T3 = zJ/3, 
the system separates into the following two phases: a low-concentration phase that is racemic 
and a high-concentration phase that is chiral and has a finite <j). High- and low-concentration 
phases correspond to an aggregate of chiral crystals and a dilute solution, respectively. 

In the following kinetic Monte Carlo simulation, we consider the case with a vanishing bond 
energy between different enantiomers, R and S; Jrs = 0. It corresponds to the interaction 
parameters K = J = Jrr. The phase diagram in the mean-field approximation turns out to 
be essentially the same as the one at K = 0, as is shown in Fig. 1(b). The tricritical point 
is now T3 = 3zJ/5. In the following simulation, the concentration of chiral molecules is fixed 
to be as low as c = 0.1. The phase separation in the mean-field approximation takes place 
at T/zJ = 0.3066 in this approximation, as is seen in Fig. 1(b). However, in low-dimensional 
systems, and for a finite-sized system in particular, thermal fluctuation is, in general, large, 
suppressing ordering. Therefore, in the following kinetic Monte Carlo simulation, a sufficiently 
low temperature is chosen to ensure that the system phase-separates with a broken crystal 
chirality. 

3. Kinetic Monte Carlo Simulation 

In a previous section, it is shown that the chiral symmetry breaking takes place as an 
equilibrium phase transition for a conglomerate system. The equilibrium argument, however, 
cannot tell us about the dynamic process, for example, how long we have to wait for the 
equilibrium state to be achieved. In the discussion on the equilibrium state, the chirality 
and position of a molecule are assumed to change instantaneously. In the actual experiment 
of organic molecules, only monomers can change their chirality in a solution. 11 Molecules 
incorporated into a cluster have to be dissolved again in a solution before they can change 
their chirality. Thus, it takes a long time to achieve homochirality. For example, in the final 
equilibrium state with homochirality at T = 0, all molecules coagulate into a single large 
crystal with a single chirality. To reach this final state, crystals have to compete with each other 
and coarsen. This coarsening process, usually called Ostwald ripening, 20 proceeds very slowly 
and thus makes the realization of chirality selection practically impossible in an observable 
time. To circumvent this obstacle, grinding plays the essential role. 10 ' 11 A dynamic simulation 
model with grinding was proposed by us previously. 18 We report here the effect of grinding 
on the temporal evolution of chirality in detail. As a measure for characterizing the dynamics 
of chiral symmetry breaking, we adopt a saturation time t s for a system that starts from a 
racemic initial state to reach the final homochiral state. We study its dependence on the total 
number of molecules N. 

Dynamics of our lattice gas model consists of three processes: (1) crystal growth and 
dissolution, (2) chirality conversion, and (3) periodic grinding. 
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3.1 Crystal growth and dissolution 

When a chiral molecule is isolated in a square lattice, it can hop to an arbitrary empty 
site at a rate k. This long-distance jump reflects the fact that the solution is strongly stirred 
and diffusion to neighboring sites has negligible effect in experiments. When the molecule 
lands on a site with some neighbors of the same chirality, bonds are formed. When there are 
n nearest neighbors of the same enatiomeric type, n homobonds are formed with the local 
energy gain —nJ. These bonds have to be broken for the next jump, and accordingly, the 
jump rate decreases to 



Since the highly coordinated molecule dissolves less frequently, homoclusters remain stable 
and grow into large crystals. In the simulation, highly coordinated molecules with n =3 and 
4 are not dissolved at all. The rate exp(— J/T) is chosen to be 10~ 3 , which corresponds to 
T/J = (31nl0) -1 = 0.1448. For the simulation at a concentration c = 0.1, this temperature 
T/zJ = 0.1448/4 = 0.036 is very low compared with the expected phase-separation line 
T/zJ ~ 0.31 in the mean-field calculation, shown in the phase diagram in Fig.l (b). At this 
temperature, the system is expected to separate into a very dilute solution and an almost 
homochiral crystal (with c ~ 1) in the mean-field approximation. 
3.2 Chirality conversion 

Experimentally, it is known that a monomer changes chirality in the solution, if an ap- 
propriate agent is added. 11 This spontaneous conversion of molecular chirality is an essential 
process for chirality selection, and we denote the conversion rate as vq\ 



when the monomers are isolated. In the following simulations, time is measured in units of 
vq and the jump rate of an isolated monomer is assumed to be as large as k = 10 4 fo, since 
the solution is vigorously stirred with a large jump rate. In a previous study, we assumed 
an enhancement of the chirality conversion rate to u n , when a monomer is in contact with n 
opposite enantiomers. The present study reveals that such an enhancement turns out to be 
inessential, as is discussed in the next section. 
3.3 Periodic grinding 

In aforementioned experiments, 10,11 many glass balls are put in the crystal growth cell, and 
when magnetic stirrers rotate to stir the solution, glass balls rotate and collide randomly with 
growing crystallites and smash them into small fragments. We simulate this grinding process 
by the following procedures:® First, the whole system is shifted randomly in space and then 
divided into square mesh cells of side length t. Each mesh is cross-cut by two diagonals into 
four triangles, and up-down or left-right triangles are exchanged randomly (see Fig. 2(a)). If 
diagonal lines cut through faceted crystal clusters, they break crystal facets and create small 
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(a) (b) 

Fig. 2. Grinding consists of (a) diagonal cut followed by up-down or left-right triangle exchange in 
each square mesh and (b) random shuffle of all meshes. 

fragments covered by high-index faces. Then, all meshes are interchanged randomly (Fig. 
2(b)). Clusters with a linear dimension larger than about 1/2 would be broken, although it is 
still possible that some clusters in neighboring meshes may reconnect to form large ones after 
grinding. Our interest is focused on the effect of grinding on the development of homochirality 
as the system size L is varied. 

The grinding process is introduced periodically in the simulation with a frequency / or 
after every interval of time, / . The relation between the grinding frequency / and the 
chirality selection is also an interesting aspect. 

4. Simulation Results 

The most interesting point we want to understand is the temporal development of chiral 
symmetry breaking. Therefore, in an initial state, equal numbers of R and S enantiomers 
are distributed randomly in the system; Nr = N$ = N/2, i.e., the initial system is racemic, 
4>(t = 0) = 0. During a simulation run, one enantiomer starts to dominate by chance and 4> 
takes a nonzero value. The system eventually reaches a homochiral state with \<f>\ = 1 at a 
saturation time t s . We expect that this saturation time t s will characterize the dynamics of 
the chirality selection. The time t s depends on many parameters, but our focus here is its 
system size dependence. For its study, the linear size of the system, L, is varied, keeping the 
chiral molecule concentration c = N/L 2 at a fixed value c = 0.1. 
4-1 Ostwald ripening 

When the simulation is performed with neither grinding nor an enhancement of chirality 
conversion rate, molecules form two large clusters of different enantiomers rather rapidly in 
our small system, and then, the competition begins between them. For example, in one MC 
simulation run of a system with a size L = 160, there are four clusters, two R's and two S's, 
at a time tuQ ~ 2.5 x 10 4 (Fig. 3(a)), and the cluster number is reduced to two with one R 
and one S each at tuo ~ 8 x 10 (Fig. 3(b)). Then, the competition of chirality conversion 



8/20 



J. Phys. Soc. Jpn. Full Paper 

starts. One enantiomer dominates slightly by chance and eats up the opposite very slowly. 
Eventually, only a single crystal cluster remains after a saturation time t s . In the example, 
tgVQ « 1.75 x 10 6 , as shown in Fig. 3(c). 




(a) (b) (c) 

Fig. 3. (Color online) Ostwald ripening: There are (a) four clusters at tv = 25,000, (b) two clusters 
at tvQ = 80, 000, and (c) a single cluster at t s VQ = 1, 750, 000. The system size is L 2 = 160 2 at a 
number of molecules N = 2560 at a concentration c = 0.1. 



What really matters is the fact that the time t s depends on the system size, as shown 
by open circles in Fig. 4. The size dependence is approximately proportional to the number 
N of chiral monomers in the system. In Fig. 4, the line t s VQ = 50CW is drawn as a guide to 
the eyes. The size dependence of t s indicates that for a macroscopic system, the homochiral 
state cannot be achieved in a practical time, even though the equilibrium argument assures 
the homochirality of the conglomerate system. The system is nonergodic. 

Here, we would like to present a qualitative argument about the time development of 
chiral order parameter eft after each enantiomer forms only a large single cluster of its kind, 
as in Fig. 3(b). Let there be Nr molecules of R-type enantiomers and N$ molecules of S-type 
enantiomers. Then, linear dimensions of clusters are approximately £r ~ \/Nr and is ~ t/Ns- 
The time evolution of the cluster size may be written as 

^5 = iff a m - - 

dt \ h £r 



dt 

with a certain appropriate kinetic coefficient K since the growth is controlled by surface 
kinetics. 22 The driving force A// is related to the concentration of enantiomeric monomers 
in the solution. Since the Ostwald ripening takes place very slowly, enantiomers change their 
chirality frequently enough so that one can assume that Afj, is the same for both chiral species, 
R and S. a represents the effect of surface tension. From the above eq. (4.1), one can derive 
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10 1 10 2 10 3 10 4 

N 

Fig. 4. (Color online) Log-log plot of the saturation time t s versus total molecular number N (a) 
without grinding (open circle), (b) with grinding but without enhancement of chirality conversion 
rate (filled box), and (c) with both grinding and rate enhancement (open up-triangle). 



the evolution equation for the size difference 

d J^l = K„(r s i-r R >). (4.2) 

Since l 2 R = Nr = N{l + 4>)/2 and i 2 s = N$ = N(l — <f>)/2, the chiral order parameter (f> evolves 
as 

# _ AKa 1 - y/T^ 

dt N <j) { ' 1 

It is solved analytically as 

y/l - (f) 2 + ln(l -y/l- (/> 2 ) = {AKa/N)t + C (4.4) 

with an integration constant C to be determined by the intial condition. Two examples of 
temporal evolution of the chiral order parameter <f> are shown in Fig. 5. The system has a size 
L 2 = 160 2 with a concentration c = 0.1, and an elapsed time is normalized by the saturation 
time t s . The saturation time t s differs by a factor of three for two runs, depending on the 
initial incubation time. However, simulation results are well fitted to the theoretical behavior 
in eq. (4.4), drawn by smooth curves, with appropriately chosen values of constant C. 
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Fig. 5. (Color online) Temporal evolution of 4> for Ostwald ripening. Two simulation results for a 
system with a size L 2 = 160 2 and N = 2560 are compared with theoretical curves, eq. (4.4), with 
an appropriate choice of C . 



By assuming a small initial cee 
2K<j(j)/N and increases exponentially as 

<f>(t) = (f)(0) exp(2K at /N) 



<C 1, the order parameter develops as dep/dt = 

(4.5) 



It then reaches the saturation value \4>(t. 

t, = 



1 at a time given by 



N 
2Ka 



In 



(4.6) 



indicating that the saturation time t s is proportional to the cluster size N, as observed in Fig. 4. 
By assuming an initial fluctuation of the order |<^(0)| ~ y/T/N \ t s may have a weak logarithmic 
size dependence; t s ~ N In N/4Ka, although it is difficult to confirm this dependence in 
numerical simulations. 
4-2 With grinding 

We introduce a grinding procedure with the frequency f /vq = 1. Then, clusters are broken 
into pieces after every time interval 1/f and competition among small clusters repeats: All 
throughout the simulation, the system contains many small clusters, as shown in Fig. 6. In a 
racemic state, there are as many R clusters as S clusters, as shown in Fig. 6(a). In a chiral 
state, the numbers of R and S clusters differ, though cluster sizes are distributed similarly. In 
the homochiral state shown in Fig. 6(b), all the clusters belong to one enantiomer, but the 
size distribution appears the same as that in Fig. 6(a). 

The chiral order parameter <j> increases steadily, as shown in Fig. 7(a). It reaches homochi- 
rality \cj)\ = 1 rapidly in comparison with the case without grinding, though the saturation 
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(a) 



(b) 



Fig. 6. (Color online) Cluster configuration with grinding (a) at the intermediate time tv$ = 1,800, 
and (b) at the saturation time of homochirality t s i/Q = 7, 200. The system size is L 2 — 160 2 at a 
number of molecules N = 2560 at a concentration c = 0.1. 
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(a) (b) (c) 

Fig. 7. (Color online) (a) Temporal evolution of cee, \<fi\, for systems of three sizes: L = 80, 120 and 
160. There are three samples for each system size. No systematic size dependence is observed for 
the initial nucleation period, (b) \<p\ close to the saturation, plotted versus t s — t, and (c) \(f>\ in a 
semilogarithmic plot. The exponential behavior exp(— (t s — i)/2500) is drawn in (b) and (c) as a 
guide to the eyes. 



time t s is distributed (Fig. 7(a)). t s does not show any systematic dependence on the system 
size L or the total molecular number N, as shown by filled boxes in Fig. 4. For a large system, 
the saturation time t s is enormously shorter than that in the case without grinding, and the 
homochirality should be established rapidly. We think that this explains essentially the exper- 
imental achievement of homochirality in a short time, observed by Viedma 10 and Noorduin 
and coworkers. 11,12 

The scattering of the saturation time t s is due to large fluctuations of incubation time 
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before the chiral order parameter <f> deviates from the initial racemic state appreciably. Incu- 
bation time is determined stochastically, and thus, the initial start-up time for cee \<p\ to shoot 
up depends on the simulation setup, more precisely stated, on the random number sequence, 
as shown in Fig. 7(a). Once cee takes a certain finite value beyond the thermal fluctuation, it 
approaches the final homochiral state in a universal manner, as shown in Fig. 7(b). In fact, 
the semilogarithmic plot in Fig. 7(c) indicates that |0| increases exponentially with time. We 
want to understand this temporal behavior. 

Without grinding, homochirality is achieved by a single large cluster that wins the com- 
petitive process of Ostwald ripening. Under grinding with a finite mesh size £, as described in 
the previous section, large clusters are broken into pieces and many clusters survive whose size 
s is distributed in a limited range, as evident in Fig. 6. From the number Cr^s(s) of clusters 
of size s of each enantiomer, one obtains the average cluster size as 

sr,s - "F^ — ~ /-r, (4-7) 

where the denominator represents the total number of enantiomeric molecules 

N r>S = J2sCr,s(s)- (4-8) 

s 

We define the total number of clusters of each enantiomeric type by 

Cr, s = Y; C R,s(s). (4.9) 

s 

An example of the total cluster size distribution C(s) = Cr(s) + Cs(s) is shown in Fig. 
8(a). It is obtained in a simulation run of a system with a size L 2 = 80 2 containing N = 640 
molecules at a time tvo = 3000, averaged over an interval of time Atis^ = 100. The grinding 
mesh size is set at I = 20. Just after the grinding, many small clusters are produced as shown 
by a curve C a (s) in Fig. 8(a), and the supersaturation is high. Then, large clusters grow 
by incorporating monomers. As time proceeds, supersaturation decreases and small clusters 
dissolve further to feed large clusters. As observed in the size distribution C&(s) after one 
period of interval 1//, just before the next grinding, the number of monomers decreases and 
small clusters disappear, whereas the number of large clusters increases, as shown by a box 
graph in Fig. 8(a). The inset of Fig. 8(a) shows that the difference AC(s) = Cb{s) — C a (s) is 
mostly positive for large sizes, indicating that large clusters grow to coarsen in this period. 
There remains exceptionally many clusters with sizes 4 and 6, since they are compact in shape 
and energetically stable. 

At the time tv$ = 3000 shown in Fig. 8(a), there is not much difference in average cluster 
size between the R and S enantiomers, sr ~ 18.9 and ss ~ 16.6, as well as in respective total 
cluster number, Cr ~ 30.2 and Cs ~ 27.6. In later time, they differ, as shown in the temporal 
evolution in Fig. 8(b). In this figure, the average cluster sizes sr and ss are shown by curves 
and the cluster populations Cr and Cs are shown by symbols. The time evolution of the 
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Fig. 8. (Color online) (a) Cluster size distribution under grinding without enhancement of chirality 
conversion rate on the cluster. Boxes represent distribution just before grinding Cb(s), and lines 
just after grinding C Q (s). Inset shows the difference between the two AC(s) = Cb(s) — C a (s), 
which is negative for small s (not shown) and positive for large s. The system size is L 2 = 80 2 
at a number of molecules N = 640 at a concentration c = 0.1. (b) Temporal evolution of average 
cluster sizes sr, s$ and total cluster numbers Cr and Cs for two enantiomeric types. The chiral 
order parameter <fi is plotted as well. 



chiral order parameter <j) is also depicted in Fig. 8(b), and the homochiral state is achieved at 
t s VQ = 8834. It is evident that as \(j)\ increases, cluster size and population start to differ for 
two species. The average cluster size, however, does not vary markedly until the very end of 
chirality selection. The size is in any case distributed rather wide, as is apparent in Fig. 8(a). 
The average size is mainly controlled by the mesh scale t and grinding frequency /. It is found 
to be independent of the system size I? . On the contrary, the cluster population increases 
proportionally to L 2 . Also, it starts to differ systematically in an early stage, even though the 
total cluster population C = Cr + Cs remains stationary, around 58 in Fig. 8(b), until shortly 
before t s . This indicates that the numbers of clusters of two enantiomers compete during the 
chirality selection. We therefore consider the selection process in terms of cluster populations, 
Cr and Cs, regarding a constant average size, sr = ss = s. 

Grinding splits large clusters to small ones, but it does not change the chirality of 
molecules. Chirality variation occurs only during the coarsening period between two suc- 
cessive grinding procedures. Therefore, let us consider the evolution of the size distribution 
{Cfl(s,t), Cs(s,t); s = 1, 2, •••} after the ith grinding until the next one, or at an interval 
of time, i/f<t<(i + l)/f. Evolution is generally described in terms of the Becker-Doering 
equation 21,22 as 

dCR d [ s,t) = -> s ) c «( s '' *) - w ( s -> s ')^( s ' *)] 

s' 
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+ 8 atl uo[Cs(l,t) - C R (l,t)], (4.10) 

where a term with Kronecker's delta 5 S) i represents the chirality conversion process for the 
monomer (s = 1). W generally represents the transition probability in size space. Then, the 
total number of R enantiomers varies as 



dN R (t) _ a v 



sC R (s,t)) 



V v(s)C R (s, t) + v [C s (l, t) - C R (1, t)] (4.11) 



s 



with the "velocity" in size space 



v (s) =J2(s' - s)W(s^ s'). (4.12) 
s' 

Since the total number of molecules N = Nr + Ns is conserved as is expressed as 

^ = [Cr(s, t) + C s (s, t)} = 0, (4.13) 

s 

v(s) should change sign as s varies. In fact, there is a critical size s c , and v(s) is negative for 
smaller sizes, s < s c , and positive for s > s c , in accordance with the sign change of AC(s) 
shown in the inset of Fig. 8(a). Those clusters with negative v(s), in fact, disappear during 
this coarsening period, as shown in Fig. 8(a). Schematically, velocity may be represented as 

v(s) = Kyfi(AfM - as~ 1/2 ), (4.14) 

with s c = (a I Afi) 2 . The supersaturation Afi, in fact, depends on the instantaneous monomer 
density. 

For the difference between the total numbers of two enantiomers, M = Nr — Ns, we have 

^=J2v(s)[C R (s,t)-C s (s,t)] 

s 

+ 2vo[C s (l,t)-C R (l,t)]. (4.15) 

The racemic state specified by 

N 

N R = N S = -, C R (s, t) = C s (s, t) (4.16) 

is a stationary solution of eq.(4.15). If the total numbers of two enantiomers differ, such that 
M = Nr — N$ > for instance, how will this difference evolve? At the end of the coarsening 
period, each enantiomer shows a cluster size distribution similar to that shown by a box graph 
in Fig. 8(a). The number of small clusters s < s c with large negative velocities v(s) vanish 
except the number of monomers. Since monomers changed their chirality often enough, their 
numbers turn out to be close for both enantiomers: Cr(1) ~ Cs(l). The numbers of clusters 
of sizes 4 and 6 remain large, which is due to the smallness of their velocities in size space. 
Their contribution to the summation on the right-hand side of eq.(4.15) is negligible. Thus, the 
amplification rate of number difference dM/dt is governed by the contribution from the cluster 
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size region with a positive velocity v(s). Then, the rate dM/dt is positive if Cr(s) > Cs(s) 
for s > s c . By assuming that the dominant contribution originates from those terms around 
the average cluster size s, the velocity v(s) is replaced by v(s) > 0, and then, the right-hand 
side of eq.(4.15) is transformed as 

^ = v(s)(C R - C s )\ t=ii+1)/f = ^V+D, (4.17) 

where Cr and C5 are the total numbers of clusters, defined in eq.(4.9) and evaluated just 
before the (i + l)th grinding. In the last equality, the approximation sr = s~s = s is used, and 
jVf(H-i) = s(Cr — Cs)\t=(i+i)/ f is the difference in number between the R and S enantiomeric 
molecules before the (i + l)th grinding. 

By assuming that the rate dM/dt remains to be of the order given by eq.(4.17) in almost 
the entire coarsening period between the consecutive ith and (i + l)th grinding processes, then 
by integration over a period 1//, one obtains 

M (*'+i) _ M (i) = ^£) M (i+i). (4.1 8 ) 
fs 

Note that grinding does not affect the total numbers of enantiomeric molecules Nr and N$, 
and thus, the value of M after the ith grinding is the same as that before the grinding M®. 
Since the chiral order parameter is written in terms of the number difference M as (ft = M/N, 
it increases exponentially after ith grinding at the time t = ij f as 

v{s)' 
sf 

with 



jW = (l - ^) V 0) = <^ (0) exp(t/r) (4.19) 



flog(l-m) «(«)' 



(4.20) 

s f . 

Thus, the homochirality \<ft\ = 1 should be accomplished at t s = — rln^A ), irrelevant of the 
system size, in agreement with the results shown in Fig. 4. 

On simply decreasing the grinding frequency from / = lto/ = 0.1, the average cluster 
size increases from s = 18, as shown in Fig. 8, to s ~ 24 and the saturation time varies 
from t s vo = 8833 to t s VQ = 15173: The larger the average cluster size s, the longer the 
saturation time t s , not in contradiction to eq. (4.20). However, the precise behavior of t s 
requires knowledge on v(s) and other factors, which have been difficult to evaluate thus far. 
4-3 With grinding and an enhancement of chirality conversion rate 

In addition to grinding, we can incorporate the enhancement of chirality conversion rate 
when an isolated enantiomer is in contact with clusters of opposite chirality. For example, the 
conversion rate is enhanced to v n = 10™fo when a monomer is in contact with n opposite 
enantiomers. n = 1 indicates that the monomer is on the flat face of the other cluster, and 
re = 2 at the kink site. The enhancement in chirality conversion actually speeds up the 
establishment of homochirality, as is depicted in the plot of saturation time t s by up-triangles 
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in Fig. 4. The saturation time t s is by far shorter than that in the case without enhancement. 
If, for instance, we include such enhancement in our sample system with L 2 = 80 2 at vq/ f = 1, 
the saturation time shortens to t s VQ = 590, instead of 8834 without enhancement. 

In our previous study, we have simulated only for a short time and thus concluded that the 
enhancement is essential for achieving homochirality. Results of the present analysis, however, 
clearly indicates that the enhancement speeds up the process but does not affect essential 
features qualitatively. The grinding alone is sufficient for establishing homochirality within 
a saturation time, which is independent of the system size. In fact, when the variations in 
average cluster size SR t s and population Cr^s are plotted with the scaled time t/t s , a universal 
feature reveals itself. Symbols of average sizes SR t g (Fig. 9(a)) and cluster populations C^s 
(Fig. 9b) collapse on those obtained without enhancement, shown by curves, in the previous 
section. 




(a) (b) 

Fig. 9. (Color online) Scaling plot of temporal evolution of (a) average sizes sr : s and (b) cluster 
populations Cr^s of R and S clusters under grinding with (symbols) and without (curves) en- 
hancement of chirality conversion rate on the cluster of opposite chirality. The system size is 
L 2 = 80 2 at a number of molecules N = 640 at a concentration c = 0.1. The grinding frequency 
is vo/f = 1. 



5. Conclusion 

A lattice gas model for conglomerate growth is shown to be reduced to a spin-one Ising 
model, the so-called Blume-Emery-Griffiths model. 19 Then, chiral symmetry breaking is un- 
destood as an equilibrium phase transition. However, the establishment of the true equilibrium 
state with homochirality is hampered by the very slow process of Ostwald ripening: Two large 
clusters of different enantiomers compete with each other, but the driving force for chiral- 
ity selection reduces its strength as the system size increases. In our simple two-dimensional 
model where crystal growth is controlled by surface kinetics, the time necessary to achieve 
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homochirality is proportional to the total number of molecules. In a macroscopic system, one 
cannot effectively reach the homochiral state with a single large crystal cluster. 12 If diffusion 
comes into play instead of stirring, the necessary time should be longer. In actual experiments 
on solution growth, the system usually ends up with many crystal clusters with different 
chiralities. 8 ' 9 ' 12 

The grinding process forcibly breaks crystal clusters, and cluster size is distributed in a 
certain range. The width of the range and the average size of clusters are almost independent 
of the system size, since they are essentially determined by the grinding frequency and mesh 
size. To compensate for the size restriction, the number of clusters becomes proportional 
to the system size. Under perpetual grinding, monomers and small clusters are constantly 
supplied by recycling and high supersaturation is maintained, since the system is closed. It 
allows the growth rate of large clusters to stay positive. The positive growth rate in size space 
amplifies the difference in total cluster number between two enantiomers, and the difference 
increases exponentially with time. The saturation time required to establish the homochiral 
state depends on the average cluster size, but it is independent of the system size. The result 
is consistent with that of experiments with grinding where homochiral state is achieved within 
an observable time. 10-12 

The enhancement of chirality conversion of monomers at the cluster surface shortens the 
saturation time, but the essential feature of chirality selection does not require enhancement, 
contrary to the conclusion in our previous paper. 18 

In a recent paper by Tsogoeva and coworkers, 14 conglomerate crystallization of chiral or- 
ganic molecule is studied. The molecule changes its chirality through a Mannich type reversible 
reaction. By assuming that the reaction proceeds fast enough so that a local equilibrium be- 
tween achiral reactants and chiral products is established, we can eliminate the degrees of 
freedom of reactants. Then, chirality conversion reduces to the reaction process Eq.(3.2), and 
our analysis possibly applies to this case. For general cases, further study is necessary. 

Grinding the system means that one has to put energy into the system from outside to 
break bondings between molecules. The system is closed but driven energetically. Therefore, 
even though the conglomerating interaction assures chiral symmetry breaking in thermal 
equilibrium, the eventual homochiral state is far away from the equilibrium state: Final steady 
state consists of many small clusters, instead of a single large crystal corresponding to the 
thermal equilibrium state. Is there any relevance of the symmetry breaking in crystal chirality 
to the homochirality in life, notwithstanding that the equilibrium state of the chiral crystal 
growth is homochiral at low temperatures while that of the chemical reaction system in life 
may be racemic? The present work shows that the homochiral state in crystal growth is 
actually a nonequilibrium state achieved and sustained by external driving. This external 
driving together with the closedness of the system enhances the nonlinear production due 
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to supersaturation and recycling. This scenario corresponds to that proposed in ref. 6 in 
a chemical reaction system. Thus, to attain homochirality in a chemical system, it seems 
necessary to place it under external driving, or in other words, homochirality in life may be 
possible in a dissipative state. 
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